Topological superfluid in a fermionic bilayer optical lattice 
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FIG. 1: The illustration of bilayer honeycomb optical lattice. 
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In this paper, a topological superfluid phase with Chern number C — ±1, possessing gapless 
edge states and non-Abelian anyons is designed in a C — ±1 topological insulator proximity to an 
s-wave superfluid on an optical lattice with the effective gauge field and layer-dependent Zeeman 
field coupled to ultracold fermionic atoms pseudo spin. We also study its topological properties 
and calculate the phase stiffness by using the random-phase-approximation approach. Finally we 
derive the temperature of the Kosterlitz-Thouless transition by means of renormalized group theory. 
Owning to the existence of non-Abelian anyons, this C = ±1 topological superfluid may be a possible 
candidate for topological quantum computation. 



I. INTRODUCTION 

Topological quantum computation, based on the ma- 
nipulation of non-Abelian anyons [H is considered as 
an effective method to deal with decoherence in real- 
izing quantum computation. The first proposed candi- 
date is the fractional quantum Hall state at filling factor 
v = 5/2 in ultra high-mobility samples 0,3- Other pro- 
posals are based on two dimensional (2D) chiral p x + ip y 
superconductors with C = ±l(the Chern- number) topo- 
logical invariable^ and then the s-wave-superconductor- 
topological- insulator systems @, 0- Among these ap- 
proaches, accurate manipulations of single quasi-particlcs 
remains a major difficulty and new techniques arc to be 
expected to overcome this drawback. 

On the other hand, cold atoms in optical lattices are 
an extensively developing research field @, Q. Because 
one can precisely controls over the system parameters 
and defect-free properties, ultracold atoms in optical lat- 
tices provide an ideal platform to study many interesting 
physics in condensed matters [l(J El- Some pioneering 
works revealed the promising potential of applying ul- 
tracold atoms to make quantum computer and quantum 
simulator. Recently, experimental realizations of quan- 
tum many-body systems in optical lattices have led to 
a chance to simulate strongly correlated electronic sys- 
tems. By changing the intensity, phase and polarization 
of incident laser beams, one can tune the Hamiltonian pa- 
rameters including the dimension, the hopping strength 
and the particle interaction at will. 

In this paper, we propose a scenario in which a topo- 
logical phase, possessing gapless edge states and non- 
Abelian anyons, is realized by proximity effect between 
a C = ±1 topological insulator and an s-wave superfluid 
(SF) of ultracold fermionic atoms in an bilayer optical lat- 
tice with an effective gauge field and a layer- dependent 
Zeeman field generated by laser-field |F2h18| . At the be- 
ginning, we give an effective design of the bilayer Haldane 
model. Then we put two-component (two pseudo-spins) 
interacting fermions on this bilayer optical lattice with 
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fixed particle concentration. For layer- 1, the Haldane 
model of two-component fermions at 1/4 filling under a 
strong Zeeman field becomes a C = ±1 topological in- 
sulator. While for layer-2, there is no Zeeman fields, we 
get an s-wave SF state by tuning the interaction between 
fcrmion via Feshbach resonance technique. Thus due to 
the proximity effect, we get an effective C = ±1 topo- 
logical SF state, of which the 7r-nux obeys non-Aelian 
statistics and becomes a non-Aelian anyon. Thus, this 
C = ±1 topological superfluid (TSF) may be a possible 
candidate for topological quantum computation. 



The paper is organized as follows. In Sec. II, we start 
with the Hamiltonian of the interacting spinful Haldane 
model on bilayer optical lattice. In Sec. Ill, we calcu- 
late the SF order parameter with mean field approach 
and get a global phase diagram at zero temperature. In 
Sec. IV, we point out that there exists a C — ±1 TSF 
due to the proximity effect of the SF order in layer-2 on 
C = ±1 topological insulator in layer-1. In Sec. V, wc 
discuss the quantum properties of the C — ±1 TSF, in- 
cluding the statistics of the 7r-flux and the edge states. In 
Sec. VI, by using random-phase-approximation (RPA), 
we calculate the phase stiffness of the C = ±1 topologi- 
cal SF. In Sec.VII, we get the Kosterlitz-Thouless (KT) 
transition temperature by the renormalized group (RG) 
theory. Finally, the conclusions are given in Sec. VIII. 
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II. THE SPINFUL HALDANE MODEL ON A 
BILAYER OPTICAL LATTICE 

In the first step, we design a bilayer optical lattice of 
the Haldane model. In Rcfs.pjj], the monolayer optical 
lattice of the Haldane model had been proposed in the 
cold atoms with three blue detuned standing-wave lasers, 
of which the optical potential is given by 



V(x,y) 



E 

J=l,2,3 



V sin 2 [&/, (x cos 6j + y sin 6j ) + n/2] ( 1 ) 



where Vq is the potential amplitude, 9± = 7r/3, #2 = 2-7r/3, 
#3 = 0, and A;l is the optical wave vector in XY plane[19|. 
On the other hand, to design a bilayer optical lattice, the 
optical potential from the interference of two opposite- 
traveling standing- wave laser beams along the z direction 
is added as 



V(z) = V L sm 2 (k z L z) - V s sin 2 (2k z L z) 



(2) 



where Vl and Vs are the amplitudes of the long and short 
laser along the z-direction. fc£ is the optical wave vector 
in z-dircction. Thus the total optical potential of the 
bilayer honeycomb lattice in our case can be written as 



V(x,y,z) = V(x,y) + V(z). 



(3) 



See the illustration in Fig.l. Since the potential barrier 
of the optical lattice along the z direction is a double well 
(See Fig. 2), the vertical tunneling between different bi- 
layer systems is suppressed seriously, each bilayer can be 
regarded as an independent 2D honeycomb lattice. The 
positions with a minimum potential along z-dircction are 



27rn ± arccos 



number. 



where n is an integer 




FIG. 2: The optical lattice potential along z-direction as a 
double well potential for AEi 3> AE2. 

Let's calculate the tight-binding model of the bilayer 
optical lattice. Firstly we derive the hopping parameter 
t± between two layers. From Fig. 2, for < Vl < 4Vg 



and AE 1 = V L + 



(AVs-VlY 
16Vs 



^ Ajfc/ 2 - i6Vs OT 

(4Vs — Vl) — > 0, one can see that the optical lattice po- 
tential along z-direction can be approximately reduced 
into a double well potential around z = 0. Then we can 
expand V(z) at z = and get 

V(z) a ~ VL {klzf - (W S - V L ) (klzf . (4) 



We denote |0) + and |0)_ as the two nearly degenerate 
ground states of the double well in the right-hand and 
left-hand wells, respectively. The corresponding eigen- 
states of the Hamiltonian are odd and even states |0) 
and |0) o which are superposition of |0) ± such that |0) o = 

1/V^(|0) + -|0)_), and |0) e = 1/V2 (|0) + + |0)_) with 
eigenvalues E$ ± AEo/2, respectively. AEq is the split- 
ting of the energy levels due to the quantum tunneling 
effect. We identify AEq, i.e., the hopping parameter t±. 

According to the instanton approach |20l [21 1 . we obtain 
the instanton solution as 



Zd (t) = ± 



_1_ / 3(W S -V L ) 
k z L V 2 (16V S - V h ) 



tanh 



(t - r ) (5) 



and then get the e nerg y level splitting AEq that corre- 
sponds to tj_ as|20l.l2ll] 



t ± = AE /2 = 2 ] /^^u }± e-^ 
where the trapping frequency uij_ is 



(6) 



d 2 V{z) 
mdz 2 



VS(W S -V L )E Z , (7) 



with z m j n being the value of z when V(z) is a minimal 
value and the classical action (S c )± of instanton is 



(S c ) 



V2{W S -V L ) (Ws-Vl) 
(WVs-Vl) V Bjf ' 



(8) 



E 



is the recoiling energy of the atoms in the 
z-direction where m is the mass of atoms. We have set 

h = i. 

Secondly we calculate the nearest neighbor hopping 
t. The optical lattice potential on XY plane is V(x, y) 
as shown in Eq. (p} which forms a honeycomb lattice 
and simulate the Haldane Model. Substituting 9\ = ^, 



2tt 



#3 = in V(x, y), we may get 



V(x,y) 



y{3+COs[K( - 



£)]+cob[«(^: 



-)]+cos[/«a;]} 
(9) 



with k = 2kL- Around the site A or B in Fig. 3, we find 
V(x t y) = 3V 2 K 2 (x 2 +y 2 )/16-3V /4 : , i.e., 



V(x, y) =mujy/2~ 3Vb/4 



(10) 



with w = y/SVoEL, where E^ = k\/(2m). Here we 
may use the the Wentzel, Kramers and Brillouin (WKB) 
method that has been extended to two dimensions to 
estimate t for honeycomb lattice. From the Ref. [13], the 
semiclassical estimation of the tunneling amplitude reads 



t = .707(-§-) 3 / 4 e 

E^ 



-0.685, 



^ Ei 



(11) 
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FIG. 3: The instanton process from A to C that leads to a 
next nearest neighbor hopping, t' . 



Finally wc calculate the next nearest neighbor hopping 
tl . The situation here is much different from that of t and 
t j_ . The energy level splitting between two wells (denoted 
by A and C) shown in Fig. 3 is due to the quantum tun- 
neling process from A to B then to C. Now we obtain t' 
to be 



/ r- a K ,/i -1.370, 

t ~ s/2te~ s * = (|^) 3/4 e 



(12) 



The coefficient 2 in e~ 2Sa comes from twice of single in- 
stanton from one site to the nearest neighbor site. In this 
sense, we have 



t'/t ~ V2e~ Sc = 1.414e 



-0.685^/Vo/iv 



(13) 



which is always a small value. For example, if we set 
Vb = 23.806Ey', the hopping parameters are given to be 
t = 0.269^ and t> = 0.05t. 

Thus when two-component fcrmions are put into this 
bilayer honeycomb optical lattice, we may get an effec- 
tive bilayer Haldane model by applying the Raman laser 
beams. Taking the tight-binding limit, we can superpose 
the Bloch states to obtain eight sets of Wannier func- 
tions w a/b ^^{t -v A/B ) = y^exp[-r 2 /(2^)], where 

I = \fh/{muo) with luq = y/ZV^E^ /%, of which the recoil 
energy reads = h 2 k 2 L / {2m) and m is a single particle 
mass. Here A, B denote two-sublattice, a = 1, 2 denote 
the index of the layer and a =t, i denote (pseudo-)spin 
degree of freedom, respectively. Then the two-component 
fcrmions in 2D bilayer honeycomb optical lattice has a 
Hamiltonian as 1231 



{i,j),a,cr ((i,j)),a,a 
-*_L.y] 4,l,cA2,cr + h.C, 



e i<fe c t 



(14) 



where c\ a a (ci t0tta ) represents fermion creation (annihi- 
lation) operators at site i of layer-a (a = 1 or 2) for 
spin a (t or |). t (f) is the real nearest (next nearest) 
neighbor hopping amplitude, and tj_ is the interlayer cou- 
pling which is much smaller than t, i.e., t± <C t. 

denote the nearest neighbor and the next nearest 
neighbor links, respectively. The next nearest neighbor 
hopping term has a complex phase 4>ij = if; where 
the positive phase is set clockwise. To design a com- 
plex phase of the next nearest neighbor hopping for a 
two-component fcrmions generated by the gauge field on 
the optical lattice, we may apply a Raman laser beams 
in XY plane with spacial-dependent Rabi frequencies as 
n sm{k L x + f )e^ and Q cos(fc L x + |)e-^ (k L = §s) 
where a denotes the length between nearest neighbour 
lattice sites. Then we get a laser-field-generated effec- 
tive gauge field on this honeycomb optical lattice similar 
to that proposed for the monolayer honeycomb optical 
lattice in Ref . (l3| . 

In addition, we apply a layer-dependent Zccman field 
to polarize spin degree of freedom by detuning Raman 
lasers only on the fcrmions in layer- 1 as 



(15) 



For big enough Zeeman field h, this term eventually 
drives this fermion model at 1/4 filling in layer- 1 to a 
C = ±1 topological insulator with fixed chemical po- 
tential fj,i = —h. Such layer-dependent Zccman field 
can be realized by the a modulated laser wave along 
z-direction, of which the wave vector is k\ but has an 
additional phase shift to the laser beams that generate 

the optical lattice, A<^ 



2tt 



For 



this case, we always get a zero Zeeman field at layer-2 
at positions with a minimum potential along z-direction 

kf 



2irn 



arccos 



but a finite Zeeman 



field at layer- 1 at positions with a minimum potential 
along z-direction k z L z = 2irn + arccos 



Furthermore, we consider a strong interaction via Fesh- 
bach resonance technique [13, HH and contact interaction 



E 



-U > 

i.a=l,2 



(16) 



where U > is the on-site attractive interaction strength 
given by integrals over the Wannier function around site 
A or B that reads 



U 



i / dr \w(r - r A/B ) 



(17) 



2V3{[Hk 2 L b 2 /4T 1 + AB/(B-B )}JV EL 



where g is coupling constant in two dimension [26j, b is 
radius of hard-sphere potential, B is magnetic field, Bq 
is resonance magnetic field and AB is the width of the 
resonance, respectively. 
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layer- 1 as 
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FIG. 4: The phase diagram of fermions (filling factor n/ = 
1/4) with attractive interaction on layer-1 honeycomb optical 
lattice for the case of t' = 0.05r,. There exist three quantum 
phases: C — ±1 QAH state, s-wave SF state and Fermi liquid 
state (FL). In this paper we only consider the C = ±1 QAH 
state of h — 7t which is marked by the red dot. 



Hi = -t ci,l,<Ai,a 



(20) 



U n< ,i,tn< ,1,4 + ^-c 



When there are a finite s-wave pairing order parameters 
of layer-1 Ai = (ci,i,fCj,i,4,), we get the effective Hamil- 
tonian as 



Hi 



off 



(21) 



Finally we get an interacting two-component fermions 
system in 2D bilayer honeycomb optical lattice of the 
Haldane model with layer-dependent Zeeman field as[23l 
MM 



H = H^i + Hz 



Hi J + H r 



where 



"Ml H cJ,l, CT Ci,l, CT - ^2 Y 4,2,^,2," 



(18) 



(19) 



with fii and [ii denoting the chemical potentials in layer- 
1 and laycr-2, respectively. 



III. MEAN FIELD APPROACH AND GLOBAL 
PHASE DIAGRAM 



Considering a tiny interlayer coupling t±, we may use 
the mean field approach separately for each layer. When 
increasing the interaction strength, the fermionic system 
described in Eq. (fTB)) is unstable against supcrfluid (SF) 
orders that are described by Ai i2 for s-wave pairing order 
parameters as Aj = {h,!^,!^) and A 2 = (ci,2,4.Ci,2,f)- 
Due to the layer-dependent Zeeman field, the symmetry 
between two layers is broken (Ai ^ A 2 ). In this section, 
we give the mean field calculations. Due to tiny inter- 
layer coupling, we may do calculations of the models for 
each layer separately. 

Firstly, we consider the interacting Haldane model on 



The energy spectrums of the fermions in layer-1 are given 
by 

(E 1>k ) 1 = h + e lik , (E ltk ) 2 = h-e 1<k (22) 
(E hk ) 3 = h + e 2k , (E hk ) 4 = h-e 2)k 

where 



e 1)fe = ^ M 2 + ([/A) 2 + 7 £ + |&| 2 + 2/, 
e 2 , fe = v /^ 2 + ((7A) 2 + 7 2 + |a| 2 -2/. (23) 
The functions j k , f are 



16 1 = t V 3 + 2 cos(V3fc„) + 4 cos(3fc x /2) cos(V3k y /2), 
7 fc = -i'[4cos(3fc 2: /2)sin(\/3fc a /2) -2sin\/3\], 

f = 2 v / 7 , 2 [M 2 + (C/A) 2 ] + |e fe |V- (24) 

By minimizing the ground state energy we arrive at 
the following self-consistent equations 



U W [ 



V i±I^tanh(-^ll) (25) 



(•Ei,k)i<° 

1 + 7fc//.._ l . / (^.fc) 



2T 



E 

<0 ^ 



(-El,fc) 2 <0 1 '* 

(^M)3) 



tanh(- 



(Bi,fc) 3 <0 



2T 



(Bi, fc )4<0 



e 2,k 



2T 
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and 



n = l+ W [ - E — tanh(— ^) 



E 



(Bi,fc)i<0 



(26) 



(Bl,fc)2<0 



tanh(-- — — ) 



(Si,*) 3 <0 



2T 



E l-(7 fc 2 + l6| 2 )// tanh( _(^,) 4 . 



(Bl,fc)4<0 



2T 



where n is fermion density and is the number of primi- 
tive cells. By this approach we obtain the phase diagram 
given in Fig. 4. 

Next we consider the interacting Haldane model on 
layer- 2 as 

~ M E 4,2,aCi,2,a ~ ^ ^ «2,*f«2,i|.- (27) 

of which the chemical potential is set to be equal to that 
of layer- 1. Considering a finite s-wave pairing order pa- 
rameters of layer- 2 A2 = (ci,2,fCi t 2,i)i we get the effective 
Hamiltonian as 



H 



- f 2J ^2Ci,2,|C ii2 ,t 4-2, <rCi,2,<7 + ^--C 

i z, cr 

(28) 



The energy spectrums of the fermions in layer- 1 are given 
by 

(E 2 ,k) 1 = (E 2 ,k) 2 = e l k ,(E2,k) 3 = (E2,k) i = e 2 fc . 
Similarly, we get self-consistent equations as 



1 _ 1 v l + 7 f// 
U AN [ ^ e, t 



tanh(^) (29) 



2T 



E 



1 - 



tanhM^ 
e 2T ;J 



and 



■''^E ''™ ^) (30) 

By this approach we obtain the SF pairing order param- 
eter given in FigfS] 




FIG. 5: SF order parameter of the fermion with attractive 
interaction on layer-2 honeycomb optical lattice for the case 
of t' — 0.05t. There exist two quantum phases: s-wave SF 
state and metal state. In this paper we only consider the s- 
wave SF state of ^2 = —7t which is marked by the red point. 



Now the total effective Hamiltonian is given by 

H c fi = i?bi + ffzecman + H$F + H c , (31) 

where 

H SF = -^E A o(cf.a.4-^.a.t + 4, a ,t d la,l) ■ ( 32 ) 
i,a 

Thus we get four self-consistent equations of the fermions 
in both layers for Ai^ and \x\p. by minimizing the ground 
state energy. We fix the fermion filling factor Uf = 1/4 in 
layer- 1. A key point here is to keep the following chemical 
potential condition: 



/'l = /'2 



(33) 



which guarantees 1/4 filling factor for the fermions in 
layer- 1. To keep the above condition, we must tune the 
chemical potential in layer-2 by manipulating the in- 
teraction U to synchronize with fj,i in layer-1. 

From FigUl one can see that for U 7^ 0, there exist 
three quantum phases in layer-1: the Fermi liquid (FL), 
the topological insulator of C = ±1 TKNN number with 
quantized anomalous Hall (QAH) effect (we also call it 
C = ±1 QAH state) and s-wave superfluid (SF)(2g|. For 
free fermions, U = 0, there exists a critical point for 
Zeeman field: h c ~ 1.37t. For the case of large Zeeman 
field, h > h c , the ground state is C = ±1 topological 
insulator with chiral edge states; For the case of smaller 
Zeeman field, h < h Cl the ground state is a FL state. 
When considering the interaction term, there may exist 
s-wave SF state, and the quantum phase transition from 
C = ±1 QAH state to SF state is the first order. In 
this paper we consider the C = ±1 QAH state with h — 
It. For this case, due to big imbalanccd Fermi gaps of 
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different spin components, the attracting interaction is 
irrelevant. The ground state is C = ±1 QAH state until 
the interaction strength is larger than 

For layer-2, the situation is much different due to zero 
Zeeman field. Fig[5]shows the SF pairing order parameter 
A2 of layer-2 for the case of t' = 0.05<, (12 = —It. For 
the weak interaction case, the ground state is FL; for the 
large interaction case, when U > 13. Oi, the ground state 
becomes SF with trivial topological properties. Just for 
this reason, we cannot get a C = ±1 topological SF in 
a monolayer system, and thus have to turn to a bilayer 
system. 

Now, we consider the model with the parameters U = 
15t, h = It, t' = 0.05t (marked by the red spots in FigH 
and FigEl respectively), at which we derived Ai = 0, 
A2 = 0.227, fii = H2 = —h = —It at zero temperature. 
Thus we have a C = ±1 QAH in layer- 1 and an s-wave 
SF state in layer-2. Now the total Hamiltonian is given 
by 



H cS = -t J2 ia^a-t' eife£ U, C 



(34) 



tj_J2(cll,a d i,2,a + h.C.) + h ^ 4,l,o Cr ?,o' 
i.a i,a,(7' 

h 53 4,i,<Ai.* + ^ 53 fi U<A2,<x 
u 53 A 2 (c l , 2 ,;c l , 2it + cj j2it at )3i4 _). 



At high temperature, the SF pairing order disappears. 
The transition temperature of SF order in layer-2 is 
k B T c ~ 2.25*. 



IV. C = ±1 TOPOLOGICAL SF DUE TO 
PROXIMITY EFFECT OF SF ORDER IN 
LAYER-2 ON C = ±1 QAH IN LAYER- 1 

We use the purterbative theory to calculate the prox- 
imity effect of the SF order in layer-2 on the C = ±1 QAH 
inlayer-1. See the illustration in Fig(5] The Hamiltonian 
has a form as 



H e fi = H + Hj 



(35) 



in which 



Ho — t 53 ^i,a,cr^j,a,o- t 53 6 c\ a <y Cj^ a ^ 
{i,j),a,a ((i,j»,a,CT 

(36) 

+ /l 53 6 !l,a<a'Ci,l,<7' ~M 53 4,1,0^,1,0 
i,(T,cr' i,<7 

~U^2 A2 ) i(Ci,2 ) J.Ci 1 2,t + 4,2,t4,2,4-) + h - C - 
i 

-^534,2,0^,2,0 



is the unperturbation term, and due to t± <C U A 2 and 
the interlayer coupling 



Hi = -t± 53 (4,i )0 -Ci,2,o + h.c. 



(37) 



is the small perturbation term. 

In the purterbative theory, we firstly use the path- 
integral representation 



Z = J [dc\ a dci^][dcl a dc2,a\e 



-So-S' 



(38) 



by replacing electronic operators c] a and Cj. a to Grass- 
mann variables c* a and c JiCT . Sq is the action as 

So = I drCo (39) 



and the Lagrangian in terms of Grassmann variables c* 
and Ci j(T is then obtained as 

£o = 53 C *i,a.ad T Ci, a .a + H (c*,c), 



where Hq (c* , c) is obtained by replacing operators in Hq 
with Grassman variables. S is the action as 



S' = / drC! 



where 



(40) 



(41) 



C! = -t± 53 ( c *,l,o C i,2,o + h.C.) . 
i,<7 

Now we integrate ± a c^^.a and get 

Z = J [dcl a dc hrJ }e-^°~ s '^ (42) 



where 



(Sl)o = J ^[534,1,0^^,1,0 - * 53 C i,l,o c 3,l,o (43) 
i,a (ij),a 

- m 53 c *,i,o c i,i,o - f 53 e ' 0i3c i,i,o c j,i,o 

+ h}]cl lta .a z Ci t i t<7 + h.c] 



and 



where ( e 



S' cS = -ln(e~ s ' 



(44) 



-s' 



J[(fc 2 <ic 2 ]e ( 52 )o s ' . Thus we have 



(e~ S ') 2 = J[dcl a dc 2 , a ]e-^oS' 



(45) 
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C=l quantum anomalous Hall state 







Interlayer 
coupling 











S-wave superconducting state 
with trivial topological properties 
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FIG. 6: The illustration of C = ±1 topological SF due to 
proximity effect between SF order in layer-2 and C = ±1 
QAH in layer- 1. 



FIG. 7: The armchair edge state of C = ±1 QAH with induced 
SF order on layer-1 and related parameters U = 15t, h — 7t, 
in = -It, t' = 0.05i, Ai, inducc = 0.001. 



Due to (S") 2 = we derive 



ft.— >-<.-'), = -i<(sT), (46) 
= J dr[^2 Ai 1 induceC»,i4Ci ) i^ + h.c] 



where 



Ai i n di; 



(47) 



t/A.,,2 

To derive this result we have used the following equation, 



( S f) 2 = ( (t-L £ C M,^^ + ^c-) 2 ) (48) 

\ i.tr / 2 

2 



£[cU,4.Ci,l,t ( C i,2,t C i,2^) : 

+ c i,i,t c *,i4 (ci,24C ii2 ,t) 2 ] 

E ku^Lt + <i,t<i,J 



V2,i 



where 



(ci,2,4.c,, 2 ,t) 2 - (Ci )2j -fC* i2) 4.) s 



£/A 



i,2 



(49) 



That means although there is no SF pairing order pa- 
rameter of fermions in layer-1, the tiny interlayer hopping 
will lead to an induced SF order due to the proximity ef- 
fect. After integrating gapped fermions on layer-2, the 
low energy effective model of such bilayer system is fi- 



nally reduced into 

//:..:: / E 4l,<Al,<r ~ *' £ ^ 



c i,l,a c j,l,<7 

(50) 



i^a.a' 2,(7 
+ E Al,induceCi,l,4.Ci,l,t + ^' C - 

where Ai induce is the induced SF order parameter of 
fermions in layer-1 which is estimated by the perturba- 

(t ) 2 

tion approach as A i induce — And for a tiny 

induced SF order parameter, we always get a really large 
energy gap of the fermions in the bulk, AE ~ 0.52i, 
which protects the topological properties of the C = ±1 
topological SF order. 

Thus due to the proximity effect between C = ±1 QAH 
and s-wave SF, the ground state is really a C = ±1 topo- 
logical SF as 



C = ±1 QAH + s-wave SF -> C = ±1 TSF 



(51) 



of which the topological properties is similar to that of 
2D chiral p x + ip y wave SF[30i. |31|. 



V. TOPOLOGICAL PROPERTIES 

In this section we will study its topological properties 
by calculating the edge states and the zero modes on a 
7r-fiux (vortex with half quantized " magnetic" flux) . 

In FigJTJ we show the gapless Majorana edge modes 
of this effective model on a lattice with open boundary 
condition along y-direction (armchair edge) and periodic 
boundary condition along x-dircction. 
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FIG. 8: The particle density of the zero modes of a pair of 
7r-flux in C — ±1 QAH with induced SF order parameter on a 
36 x 36 lattice in layer-1. In the right figure, there is a branch- 
cut (red dashed line) between two 7r-fluxes (cyan plaquette) 
which changes the signs of the hopping terms on green links. 

In Fig J5] wc demonstrate the existence of a single 
zero-mode on a 7r-fiux obtained as a solution of the 
Bogoliubov-de Gennes equations [HI, [33|]. This is the 
Majorana zero energy mode and can be described by a 
real fermion field 7^ = J drluotp^ + VQijj] (7^ = 7) j32|. 
When two 7r-fluxes are fused together (taken to the same 
plaquette in the honeycomb optical lattice), the result 
contains more than one quasiparticle due to the Ising 
fusion rule [Hj], 

a x cr = I + 4>. (52) 

These results imply that the 7r-fiux is a non-Abelian 
anyon (people also call it Ising anyon) obeying non- 
Abelian statistics as that in chiral p x + \p y wave SF. The 
topological properties of the C = ±1 topological SF or- 
der are similar to those of chiral p x + ip y superconductors 
with C = ±1 topological invariable. 

VI. PHASE STIFFNESS OF C = ±1 
TOPOLOGICAL SF ORDER 

In this section, by using the random-phase- 
approximation, we calculate the phase stiffness of 
C = ±1 topological SF order which is determined by 
the phase fluctuations in layer- 2. The Hamiltonian of 
fermions in layer-2 is 

H 2 ,oS = ~t ^ £ U,Jj,2,<T ~ t' ^ ^4,2,^,2,^ 

~ V X! 4,2,a^,2,cr ~ U ^ ™»,2,t"i,2,4. + h - c - ( 53 ) 
i.a i 



Thus in path-integral representation the partition func- 
tion is given by 




by replacing electronic operators c] 2 a and Cj t 2,<? to 
Grassmann variables c* 2(T and Cj-,2,<r- The effective ac- 
tion S 2 is 



S 2 = / dr£ 2 (55) 



and the Lagrangian in terms of Grassmann variables is 
then obtained as 

£2 = X] c* a<J d T c ia ,a + H 2 , e s(4 <7 ,c 2 ,a), (56) 

/.rr 

where H 2c g(c 2 a , C2,a) is obtained by replacing the oper- 
ators in H 2c r with Grassman variables. 

Considering the s-wavc pairing order parameter of 
layer-2 A2 = (02^.-1-^2,^4), we get the partition function 
as 



Z= / [dc* 2 a dc 2 ^][dA*dA]e- s 1 (57) 



where the action of fermions in layer-2 is 
S = ( iw ™ + 7k - v)ck, a Ck,a + 2/3//JV (58) 

uj m ,k 

+ &C* A M,a c BM,a + C.C. 

~ 77^ H( f/A *( _ «) c - k +q/ 2 ^t c k+q/2.i 

+ ^A(?)4 +q /2^ c -k +q /2, t ) + 2 ^E UA *(i)Mq) 

with uj m = (2n + 1)tt/(3 and wj = 2mr//3. We denote 
k = (\ijj m ,— k) and q = (iwi,— q). After integrating over 
the fermionic field, the action turns into 

S = 2/3A^?7A*(q)A(q)+2/3/xAr- ^ Trln(G _1 ), 

q k,k'.q 

(59) 

where 
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G~ 



-UA*(-q) 




-UA(q) 

+ 7k + A^k.k' 





6A,k' 



(iw m - 7k - At)^k,k' 

-C/A*(- 9 ) (iw„ 





-Cfe^k.k' 

- 7k + M)£k,k' 



Then we consider the phase fluctuations on the SF or- where 
der parameter and define UA(q) = A + A(q) where the 
mean field value Ao is a real constant and A(q) is com- 
plex fluctuating field. Thus G -1 is divided into two parts, 
i.e., 

G~ 1 = Gq 1 +G1 1 (60) 



1 



i + 7k 

-Ao 

o 



-Ao 
iuj m + 7k + n 


-ft 



o 

- 7k 
-Ao 



\ 





-& 
-A 

- 7k + M) / 



and 








-A B (-g) 




Using the expansion of the logarithm, one gets 



TrlnG" 1 = Tr ln^ 1 + Gr 1 ) 

f 



(61) 
(62) 



Trln G n 



Tr(GoGr GqG7 ) 



Next we investigate the Gaussian fluctuations of the 
paring field A(q) around the saddle point A [H, [36| . The 
fluctuation field is written as 

A\q) = (A A (-q),A A (q),A B (q),A* B (-q)). (63) 

Then quadratic effective action becomes 

S = S + Si + 2(3N UA * (q) A (q) + 2/3/iiV, (64) 



where 



5o = ^TrlnG 1 , 



5 i= 2^^ At(q)Q(q ' k)A(q) ' 



(65) 
(66) 



q,k 



and the detailed form of elements in Q(q, k) are shown in 
Appendix A. Using the Matsubara summation formula, 



one can obtain the quantities Qij above. Then in the 
static limit, i.e., iw; — > 0, at zero temperature, Qij can 
be described as follows, for example, if 



Qi,j = 2J G Qbc (k)G gh(k - q) 



(67) 



and then after the summation, it becomes 

n A bc (k)B gh {k - q) + B bc (k)A gh ( k - q) 

lij 



ei,k- 


q + ei,k 




A bc (k)D gh {k-q) -\ 


-B bc (k)C gh (k- 


q) 


e 2,k-q 


+ ei,k 




C bc (k)B gh (k-q)l 


D bc {k)A gh (k- 


q) 


ei,k- q 


+ e 2 ,k 




G 6c (k)ZV(k-q)7 


-D bc {k)C gh (k- 


q) 


e 2,k-q 


+ e 2 ,k 





(68) 



where the parameters Aij , Bij , dj , Dij are all shown in 
Appendix A. 

In order to obtain the superfluid phase stiffness, we 
further separate the fluctuation into its amplitude and 
phase components A A/B (q) = [A A/B (q) + i9 A/B (q)]/V2 
with real filed A^/s(q) and 9 A / B (q). The changes of 
basis can be written as 



/ A^(-q) 

Ayi(q) 

A B (q) 
Ui(-q) 




/ Aa(cj) 
^A(q) 

Mq) 
V A B (q) 



(69) 
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Then we have 



q 



for q 2 0. 

After obtaining the phase stiffness of C = ±1 topo- 
(70) logical SF order, the effective Lagrangian of the phase 
fluctuations is obtained as 



where A = ( A^(q) 0,4 (q) #s(q) As(q) ) , and the de- 
tailed forms of elements in matrix W are shown in Ap- 
pendix A. Integrating over the gapped field A(q), in the 
static limit at zero temperature, we obtain 



Si[0] = §£[Mq),Mq)]* 



(71) 



= ^.(O)(V0) a 



(78) 



For example, for the case of U = 15t, t' = 0.05£, 
fi2 = —h = — It, we have a small phase stiffness as 

p s (0) ~ 0.00286i. 



with 



X 



W 22 

w 32 



Wbu W23 
W B 21 W 33 



W B 12 
W B 22 



(72) 



where the elements Wbu, Wbi2, Wb21, Wb22 become 



Wbu = W A iiW2iW 12 
+ WA21W24W12 - 

W B 12 = WAI1W21W13 

+ W A 2lW 24 W 13 - 
W B 21 = W A uW 3 lWl2 

+ Wa2iW 34 W 12 - 
W B2 2 = W A uW 31 W 13 
+ W A 2iW M W 13 - 



- WA12W21W42 
WA22W2 4 W 4 2 

- WA12W21W43 
W A 22W 24 W 43 

■ W A i2W 31 W i2 

W A 22W 34 W 4 2 

■ W A i2W 31 W i3 

W A 22W 34 W 43 



(73) 



and the elements Wau, Wai2, Wa21, Wa22 are given by 

W44 



Wau 

W A \2 
WA2\ 
WA22 



W 41 W 44 - W 14 W 14 

-w 14 

WuWa-WiiWu 

-w 41 

W n W 44 - W 14 W 14 
Wn 

W n W 44 - WuWu 



(74) 



At last, we arrive at the the effective action for phase 
fluctuations as 



S[0] = § £[Mq),MqF 
q 

where the elements of matrix T read 



0b(cl) 



Tn(q) = i 



Ti 2 (q) 
Tai (q) 



r 22 (q) = 



(W22 - Wbu) 

{W 23 - W B 12) 

(W32 - WB2X) 

(W33 - W B 22) 



(75) 



(76) 



We may derive the zero temperature superfluid stiffness 
p s (0) numerically in the static limit by identifying [35T - [37| 



>/Tii(q)T 22 (q) - ^I2(q)?2i(q) 



(77) 



VII. KOSTERLITZ-THOULESS TRANSITION 

From above calculations, one can see that the induced 
SF pairing in layer- 1 will be determined by the SF pair- 

ing in layer-2 as Ai^nducc = — ■ If there exists a 
vortex in layer-2, A 2j i — > A2 exp[i J^i^a Imm(z, — zi)], 
there will appear induced vortex in layer-1, Ai : j n d U cc,i — > 
Aia n duccexp[-i^/^ilm ]n(zi - zi)] where z l is the posi- 
tion as zi = Xi + iyi- So we can only study the dynamics 
of vortices in layer-1 which is defined as 



Imln(^ 



zi) 



(79) 



where Zi is the vortex position as Zi = x, + iyi. 

According to the above analysis, we can get the 
Kosterlitz-Thouless (KT) transition temperature by the 
renormalized group (RG) theory. For two vortices in 
layer-1, there exists a confinement potential as 



V ~ q 2 In ■ 



(80) 



at |r| >a, where q 2 = 27T/o s (0), and r is the distance 
between the vortex and anti-vortex. With the increase 
of temperature, the vortex-antivortex pairs can be ther- 
mally excited, leading to a contribution to the screening 
effect by reducing V to V g = -V, where K denotes the 
dielectric constant. In the following, we shall treat the 
screening effect based on an RG treatment. 

In an RG procedure, the contributions from the pairs 
with the sizes between r and r + dr will be integrated 
out, starting from r = a. The probability for the vortex- 
antivortex pairs separated by a distance r is controlled 
by the pair fugacity y 2 {r). In the KT theory [H, l39j . 
the initial is y 2 (a) = e~^ Ea , (where j3 = ^Kj> and E c 
is the core energy). The renormalization effect is then 
represented by two renormalized quantities, X(r) = 
and y 2 {r), which satisfy the following famous recursion 
relations 



7T 



dy/dl=(2--)y, 



dX/dl 



4ir 3 y 2 , 



(81) 
(82) 
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where r = ae. From Eqs. ([81j) - (l821) . we find 
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1 



+ —(X-X ) 



1 1 X 



(83) 



where X = X(l = 0) = -f^ (with = 0) = 1). The 
RG flow is then obtained from Eq. (f52"|) by 



dX' 



-Y, 



47r 3 F 2 + 4(X' - X ) - 2tt\u{X'/X ) ' 



(84) 



The pair fugacity can be determined by y 2 (l) = 
e -2/ '(2-fK_ 

The RG flow diagram of Eqs.(JHI])-([82]) is as follows: 
the two basins of attraction are separated by the initial 
values which flow to X* — > ^ and y* — ► in the limit 
I — >• oo. In terms of Eq. (f84| . the separatrix of the RG 
flows is given by 



Haldane model of two-component fermions at 1 /4 filling 
under a strong Zeeman field becomes a C = ±1 topologi- 
cal insulator. While for layer-2, there is no Zeeman fields, 
we get an s-wave SF state by tuning the interaction be- 
tween fermion via Feshbach resonance technique. Thus 
due to the proximity effect, we get an effective C = ±1 
TSF state. We also study its topological properties and 
then show the gaplcss Majorana edge modes and the non- 
Abelian statistics of the 7r-flux. This C — ±1 TSF there- 
fore may be a possible candidate for topological quantum 
computation. Finally we calculate the phase stiffness by 
using the RPA approach and then derive the temperature 
of the KT transition for the system. 



Acknowledgments 



dX' 



Xo 4(X'-f)-27rln(2X'/7r) 



(85) 



Based on the RG equation of ([8"5"j) . one can determine the 
KT temperature Xkt- Finally we approximately have 



X(l) ~ X 



xa -> oo) = - 

v ; 2 



and 



k B T KT — — = 



tt Ps (0) 



(86) 



(87) 



For the case of U = 15i, t' = 0.05t, \i\ = \x-z = —h = 
-It, we found a fairly low KT transition temperature as 



B ^ KT 



7Tp s (0) 



0.0045t. 



(88) 



Below Tkt — 0.0045t/fcs, we have a TSF with long range 
phase coherence. In the temperature region Tkt < T < 
T c ~ 2.25t/kB, we have the SF pairing but no phase 
coherence. And in this region, the vortex arc deconfined 
from the bound state. At higher temperature, T > T c ~ 
2.25t/ks, the SF pairing order disappears. 
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Appendix A: Parameters in Green functions 



In the appendix, we first give the elements in Q(q, k): 



Qi2 = ^Goi 2 (fc)Goi 2 (£;-<?) 

Ql3 = ^Goi4(fc)G 03 2(fc-9) 

Mm 

Ql4 = 4 X] G 013(fc)G 42(fc _ 



(Al) 



VIII. CONCLUSION 



In the end, we conclude our discussions. We propose 
a scenario in which a topological phase, possessing gap- 
less edge states and non-Abclian anyons, is realized by 
proximity effect between a C = ±1 topological insula- 
tor and an s-wave SF of ultracold fermionic atoms in a 
bilayer optical lattice with an effective gauge field and 
a layer-dependent Zeeman field generated by laser-field. 
At the beginning, we give an effective design of the bi- 
layer Haldane model. Then we put two-component (two 
pseudo-spins) interacting fermions on this bilayer optical 
lattice with fixed particle concentration. For layer- 1, the 



Qn = \Y, G ° 21 ( fc ) G ° 21 ( fc - «) (A2) 
Q 2 2 = ^G 22(A;)Goii(fc-<?) 



/23 



^G 02 4(A;)Go3i(fc-<?) 
^X)G<H3(*)GW(*-3) 
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14 



and 



with 



A 44 (k) = 

B 44 (k) = 
C 44 (k) = 

#44 (k) = 

p 44 (k) 
g 44 (k) 



1 


4- vaa (k) 


'Ylr - 

/K 


h U + 






4 




4pi 




1 




/k 




y 44 v*v 




4 




4e 4 




1 


-p 44 (k) 


-7k - 




? 44 (k) 




4 




4e 2 




1 


-p 44 (k) 


-7k - 


hM- 


g 44 (k) 




4 




4e 2 





(A30) 



-M7k 



VTkV + AfTTkklV 

- 7k (/x 2 + A§ )+^( 7 g + |g k p 
V7^[M 2 + A^] + iaiV 



(A31) 



The detailed forms of elements in matrix M are given 
as follows: 



V^i = ^^(Qii+Q 2 i+Qi2 + Q22) (A32) 

k 

W 12 = 2^ E ^11 + ^l) - + ^ 22 ) 



11/13 = 2JV ^ + Q 23) + *(Ql4 + Q2 4 ) 



^14 = ^ 2(0l3 + <923 + Ql4 + Q 24 ) 



^21 = ^ X - Oai) - z (^i2 - ^22) (A33) 

k 

2^7 X^ 11 - Q 2i - Q12 + Q22) 

k 

^2(-Qis + Q23 + Qia - Q 24 ) 

k 

W 2 4 = — -i(Qi3 - Q23) - «(Qi4 - Q 24 ) 



W 22 

1T 23 



^1 = E ^31 - Q41) + *(Q 32 - Q 42 ) (A34) 

k 

^2 = J2(-Qsx + Q32 + Q41 - Q 42 ) 

k 

W33 = — E(Q 3 3 - Q43 - Q34 - Q44) 

k 

^34 = TN X l ( Q 33 ~ Q43) + j(Q 34 - Q 44 ) 



Wax = 7^ ^2(Qsi + Q41 + Q32 + 04 2 ) 
W 42 
W43 



(A35) 



k 

- ^ ~i(Qs3 + Q43) + «(<3 3 4 + Q44) 



Wii = 2^ X^ 33 + Q43 + Q34 + Q44). 
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